Longitudinal magnetic excitations in classical spin systems 
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Using spin dynamics simulations we predict the splitting of the longitudinal spin wave peak in all 
antiferromagnets with single site anisotropy into two peaks separated by twice the energy gap at 
the Brillouin zone center. This phenomenon has yet to be observed experimentally but can be easily 
investigated through neutron scattering experiments on MnF2 and FeF2. We have also determined 
that for all classical Heisenberg models the longitudinal propagative excitations are entirely multiple 
spin-wave in nature. 



The mechanism for longitudinal excitations in high 
spin magnets with weak to moderate anisotropy, like 
MnF2, FeF2, RbMnFg, EuO, and EuS, is not completely 
understood. There are conflicting theoretical predic- 
tions J iB] and experimental results of limited resolu- 



tion ||2-|]; however the spin dynamics simulation tech- 
nique is able to analyze both the transverse and longi- 
tudinal components of the dynamic structure factor for 
simple classical Heisenberg models. This is true in both 
the hydrodynamic and critical temperature regimes and, 
unlike mode coupling theory, the accuracy of our results 
can be improved continuously through the use of more 
computer time. Indeed using high speed supercomputers 
we have already achieved considerably higher precision 
than existing experimental results. 

The above materials all have spin values {S > 2) which 
are large enough to be effectively described by the classi- 
cal limit, 5* — » oo, and bi-linear exchange interactions be- 
tween nearest, and in some cases second neighbor atoms 
on simple lattice structures. RbMnFa is an SC antifer- 
romagnet, MnF2 and FeF2 are BCC anisotropic antifer- 
romagnets with weak and moderate anisotropy respec- 
tively, and EuO and EuS are ferromagnets. The degree 
of anisotropy in EuO, EuS, and RbMnFs is negligible. 
While for EuO and EuS, and to a lesser extent MnF2 
and FeF2, the second neighbor interactions are not neg- 
ligible, a qualitative understanding of the magnetic dy- 
namics can still be obtained through a model with only 
nearest neighbor interactions. 

We performed our simulation on the isotropic SC 
Heisenberg magnet with both ferro- and antiferromag- 
netic bi-linear interactions and the anisotropic BCC 
Heisenberg antiferromagnet. The Hamiltonian for our 
model is given by 



7i= J^Sr-Sr' -Dj2(.^r? 
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where ^r is a three-dimensional classical spin of unit 
length, (rr') is a nearest neighbor pair, and D is the 
uniaxial single-site anisotropy term. We determined the 
dynamic structure factor in the [100], [110], and [111] re- 



ciprocal lattice directions. For the antiferromagnetic case 
we have made the transformation q —^ q + Q where Q is 
the Brillouin zone boundary in the [111] direction. 

We have used the spin dynamics simulation technique, 
which has been developed in previous work [HQ, to cal- 
culate the dynamic structure factor. The spin dynamics 
simulation technique involves the creation of an equi- 
librium distribution of initial configurations using the 
Monte Carlo (MC) technique which are then precessed 
through constant energy dynamics to yield the space- 
time correlation function from each configuration. These 
are averaged together and Fourier transformed to obtain 
a result for the dynamic structure factor. By includ- 
ing more initial configurations the accuracy can be im- 
proved indefinitely. We have used periodic boundary con- 
ditions and studied lattices of linear sizes of L = 12 and 
L = 24. The critical temperatures were determined using 
the fourth order cumulant crossing technique tt for the 
anisotropic systems, and Tc is already accurately known 
for the isotropic Heisenberg model M . Anisotropy value 
D — 0.0591 was used for MnF2 to match the experimen- 
tally determined [^ degree of anisotropy. 

For all these models we performed the simulation at 
temperature T = 0.5Tc, low enough to be completely 
outside the critical regime but not too low for the MC 
simulation to produce an equilibrium distribution of con- 
figurations. We also studied T = O.ST^ and T = O.QT^ 
for the isotropic ferromagnet to investigate the approach 
to the critical regime. The results from the higher tem- 
perature simulations were convoluted with a Gaussian 
resolution function of width S^^ to minimize effects due 
to the finite time cutoff. For the isotropic case we used a 
distribution of 1000 initial configurations for L = 12 and 
200 configurations for L = 24, and for the anisotropic 
case we used 6000 configurations for L = 12 and 400 
configurations for L = 24. A smaller number of config- 
urations was used for L = 24 due to limits of computer 
time. 

For the longitudinal component of the dynamic struc- 
ture factor in the ferromagnetic case we observed many 
excitation peaks, as seen in Fig. nl however a different set 



is present for the L = 12 and L = 24 lattice sizes at the 
same q value. We conjecture that these excitations are 
two-spin-wave peaks since due to finite lattice size effects 
the frequencies of the peaks will be limited to certain val- 
ues which will be different for different lattice sizes. In 
order to test this hypothesis we need to be able to predict 
the expected positions of the two-spin-wave peaks. This 
requires that we obtain an approximate estimate for the 
general form of the dispersion curve at the temperature 
at which the simulation is performed, i. e. i-u{q) at all 
points on the reciprocal lattice. 
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FIG. 1. The longitudinal component of the dynamic struc- 
ture factor for the isotropic SC ferromagnet vs. frequency. 
The predicted positions of the two-spin-wave peaks are su- 
perimposed on the graph of the dynamic structure factor. 
Note that a logarithmic scale has been used for the dynamic 
structure factor. 

The spin wave stiffness coefficient, D{T), is determined 
from the low q, low temperature limit of the dispersion 
curve where lj = D(T)q^ for the ferromagnetic case and 
Lo — D(T)q for the antiferromagnetic case. We assume 
that this holds above the low temperature limit, and the 
dispersion curve will be given by the linear spin-wave 
dispersion curve, i^(q) at T = 0, multiplied by the factor 
D{T)/D{0). At T = 0.5Tc the prediction diverges from 
the real result only slightly in the high q region where the 
prediction is at a lower cu value than the actual result. 
With increasing temperature our approximation of the 
form of the dispersion curve becomes less accurate. 

Now that we have determined that the linear spin- wave 
result multipfied by the factor D{T)/D{0) is a good ap- 
proximation to the dispersion curve at T = 0.5Tc, we 
can estimate the spin- wave frequency over all q values, 
not just those in the reciprocal lattice directions we have 
measured. All two-spin-wave creation peaks will thus be 
at frequency 



.(qi±qj) = cj(qi)+cj(qj) 



(2) 



and the spin-wave annihilation peaks will thus be at fre- 
quency 



where qi and qj are the wave vectors of the two spin 
waves which comprise the two-spin-wave excitation. As 
shown in Fig. n, for the isotropic ferromagnet we see a 
near perfect match of peaks in ^(q, lu) and the predicted 
positions of the two spin- wave annihilation peaks, clearly 
indicating that the excitations in the longitudinal com- 
ponent are due to two-spin-wave annihilation. 

For the anisotropic antiferromagnet the approximation 
of the dispersion curve from a measurement of D{T) is 
inaccurate since the actual dispersion curve does not fol- 
low the functional form of the zero temperature disper- 
sion curve. Instead we identified a set of two-spin-wave 
annihilation and creation peaks which involve spin-waves 
exclusively in the reciprocal lattice directions which we 
measured. Plotting the expected frequencies of these 
two-spin-wave peaks we see a good match between the 
longitudinal spin-waves and the expected values for both 
the annihilation and creation peaks. This result for weak 
anisotropy (MnF2) is shown in Fig. g. Note that no trace 
of a peak is seen at the location of the single spin-wave 
peak in the transverse component. We see, as expected, 
the presence of both creation and annihilation spin- wave 
peaks for the isotropic antiferromagnet as well. 
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FIG. 2. The longitudinal component of the dynamic struc- 
ture factor for -the weakly anisotropic antiferromagnet, MnF2 
vs. frequency. We compare the predicted positions of 
two-spin-wave peaks to the actual peak positions. The dot- 
ted lines are predicted positions of two-spin-wave annihila- 
tion peaks and the dashed lines are the predicted positions of 
two-spin-wave creation peaks. Note, that once again we have 
used a logarithmic scale. 

The two-spin-wave peaks where one of the single spin- 
waves of which they are made up has the lowest q are 
the most intense. This is to be expected since the lower 
q single spin- waves are more intense themselves. As the 
temperature rises the two-spin-wave peaks broaden. For 
the antiferromagnetic case, where the longitudinal and 
transverse components can not be separated, as T ap- 
proaches Tc the two-spin-wave peaks disappear into the 



tails of the single spin-wave peak and the diffusive cen- 
tral peak. For the ferromagnetic case the two-spin-wave 
peaks do not disappear entirely, the two-spin-wave peak 
corresponding to the lowest q spin-waves remains and 
the other two-spin-wave peaks broaden to disappear into 
its tail. A previous study by Chen et al. Q misidentified 
the peak in the longitudinal component as a single spin- 
wave excitation. This happened because they were only 
able to look in the [100] lattice direction for which the 
dominant two-spin-wave is at the same frequency as the 
single spin-wave peak for any given q. Since for ferro- 
magnets only spin- wave annihilation peaks are present, 
the dominant two-spin-wave process found in the [100] di- 
rection is [q,0,l]-[0,0,-l]. Since in the low q limit u oc q^, 
w([g,0,0]) w a;([g,0, 1]) - a;([0,0, -1]). This is however 
not the case for the [111] direction. In Fig. g we show the 
longitudinal and transverse components in the [100] and 
[111] directions at temperatures approaching Tc- Note 
the fact that the longitudinal and transverse spin wave 
peaks are at the same frequency in the [100] direction but 
not in the [111] direction. 
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FIG. 3. The longitudinal and transverse components of the 
dynamic structure factor vs. frequency as temperature ap- 
proaches the critical regime are shown; a) shows the [100] 
direction and b) shows the [111] direction. A resolution func- 
tion with coefficient 5^ — O.Ol/jJj was used. The data are 
unnormalized to show the relative intensities of the longitu- 
dinal and transverse spin- wave peaks. 

The existence of a set of finite two-spin-wave peaks 



is an artifact of the finite size of the lattice we use in 
our simulation. If one were to measure the longitudi- 
nal component of the dynamic structure factor for a real 
crystal, i. e. effectively an infinite lattice of magnetic 
moments, the spectrum of possible two-spin-wave peaks 
would be continuous. The longitudinal component of the 
dynamic structure factor for the isotropic antiferromag- 
net has been measured experimentally by Cox et al. [0. 
Taking measurements in the [111] direction, they found a 
diffusive central peak and a propagative spin-wave peak 
at the same frequency as the transverse single spin- wave 
peak. The longitudinal component of the dynamic struc- 
ture factor for the isotropic ferromagnet EuO has been 
found experimentally to have a propagative spin-wave 
and no diffusive central peak [g|. 

For both the ferro- and antifcrromagnetic cases, the 
longitudinal excitations at the single spin-wave peak fre- 
quency can be explained in terms of two-spin-wave peaks 
as follows. Since spin-wave intensity increases with de- 
creasing q the most intense two-spin-waves will be those 
which are comprised of one spin-wave with a very small 
q and another spin- wave with a q very close to the resul- 
tant two-spin-wave q. If we look at the limit of infinite 
lattice size, i.e. the real system, for the isotropic case 
a spin-wave with extremely small q will have a negligi- 
ble frequency. Thus the intensity of the two-spin-wave 
creation and annihilation peaks will have a maximum at 
the single spin-wave frequency. Even though only annihi- 
lation two-spin-wave peaks are present for the ferromag- 
netic case one should still see this effect in both the ferro- 
and antifcrromagnetic cases and this is what is seen in the 
experiments |3l0| . We see no evidence of a diffusive cen- 
tral spin- wave peak in the longitudinal component of the 
dynamic structure factor for the isotropic ferromagnet in 
agreement with the experimental results of Dietrich et 
al. [^, and the theoretical result of Villain |ll|] but in 
disagreement with the theoretical predictions of Vaks et 
al. [|. 

When wc apply this same reasoning to the anisotropic 
antifcrromagnet, which as shown in Fig. y also displays 
the same two-spin-wave peak behavior, we are left with 
an extremely intriguing result which can be measured ex- 
perimentally. A spin-wave at very small q will no longer 
have a negligible frequency but instead the frequency of 
the energy gap at the Brillouin zone center. If our hy- 
pothesis about the origin of excitations in the longitudi- 
nal component of the dynamic structure factor is correct 
then for an anisotropic antifcrromagnet one will observe 
an apparent splitting in the spin-wave peak in the longi- 
tudinal component. In an infinite lattice there will be a 
peak due to two-spin-wave creation which is shifted up- 
wards from the single spin- wave frequency by an amount 
equal to the energy gap, and a spin- wave annihilation 
peak shifted downwards by the same frequency. 

The theoretical explanation as to why the longitudinal 
component of the dynamic structure factor is made up of 



two-spin-wave peaks and why only two-spin- wave annihi- 
lation peaks are present for the ferromagnetic case while 
both two-spin-wave annihilation and creation peaks are 
present for the antiferromagnetic case is as follows. If we 
express the dynamics in quantum mechanical formalism 
the spins, as they appear in the Hamiltonian, are ex- 
pressed in terms of the operators S^, S^ , and Sf where 






(4) 
(5) 



S^ and S~ can be expressed in terms of ladder operators, 
Oi and a^, which raise or lower Sf by one quanta. In the 
linear approximation 



sr = (2S)i/2a^. 



(6) 

(7) 



If we take the approximation for Sf to one order higher 
than Sf = S then 



a]ai 
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(8) 



For the classical limit this corresponds to an infinitesimal 
raising and lowering of the z component of the spin which 
is a longitudinal spin-wave. These spin-waves will be at 
the frequencies corresponding to the difference between 
the frequencies of single transverse spin-waves since they 
result from creation followed by annihilation of a spin- 
wave. Thus for the ferromagnetic case all two-spin wave 
excitations will result from these annihilation processes. 

Unlike the ferromagnetic the antiferromagnet is not 
a quantum state of the Heisenberg Hamiltonian. As 
a result the a.; and aj' are replaced through a Bogoli- 
ubov transformation by new operators which arc a linear 
combination of the creation and annihilation operators 
but which only connect excitations on the same sublat- 
tice ||l^ . As a result both creation and annihilation two- 
spin- wave excitations exist. 

In conclusion, for magnets which can be described by 
a classical spin model, the longitudinal propagative ex- 
citations are made up of two-spin-wave peaks. In agree- 
ment with our theoretical picture only annihilation two- 
spin-wave peaks were present for the ferromagnetic case 
but both annihilation and creation two-spin-wave peaks 
are present for both the isotropic and anisotropic anti- 
ferromagnets. As the temperature approached the crit- 
ical regime the longitudinal two-spin-wave peaks broad- 
ened and converged into a single peak at the frequency 
of the dominant lowest-g two-spin-wave peak. In an 
isotropic lattice of infinite size, i.e. a real crystal, the 
two-spin-waves will result in a peak at the single spin- 
wave frequency for either ferromagnetic or antiferromag- 
netic isotropic systems, in agreement with experimental 
results [pP]. For the anisotropic antiferromagnet both 
creation and annihilation two-spin- wave peaks exist, and 



the longitudinal component of the anisotropic antiferro- 
magnet should show two peaks; a two-spin-wave annihi- 
lation peak at the single spin- wave peak frequency minus 
the energy gap frequency, and a two-spin-wave creation 
peak at the single spin- wave peak frequency plus the en- 
ergy gap frequency. This result indicates the presence of 
a new form of excitation behavior in magnetic materials 
which can be directly tested experimentally. 
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